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Abstract 

We calculate the drag force on a heavy quark hit by a shock wave, thus generalizing the 
strongly coupled AdS/CFT heavy quark drag calculations to both hot and cold nuclear 
matter. The derivation employs the trailing string configuration, similar to that used in 
the literature for a quark moving through a thermal medium, though in the shock metric 
the string profile is described by a much simpler analytic function. Our expression for 
the drag depends on the typical transverse momentum scale of the matter in the shock. 
For a thermal medium this scale becomes proportional to the temperature, making our 
drag coefficient and momentum limit of applicability identical to those found previously. 
As the shock wave can be composed of either thermalized or non-thermalized media, our 
derivation extends the existing drag calculations to the case of arbitrarily distributed mat- 
Q_i ter. 
<D 

i— —i Keywords: AdS/CFT correspondence, heavy quark energy loss, jet quenching 

^ PACS: 11.25.Tq, 12.38.Mh, 24.85.+p, 25.75.-q 

> 
CO 

1 Introduction 

O The simultaneous emergence of the anti de-Sitter/conformal field theory (AdS / CFT) correspon- 
dence [1-4] as a tool to study strongly coupled systems (see [5, 6] for reviews) and the failure 
of perturbative QCD (pQCD) techniques [7] to quantitatively describe the observed heavy ion 
physics phenomena has inspired both a critical reevaluation of the traditional methods and an 
explosion of new research [8-14]. In particular a detailed theoretical understanding of high- 
er particles and jets, coupled with experimental measurements associated with them, holds 
out the promise of uniquely probing the high-pr physics, low-p-r collective physics, and their 
mutual interaction in a bulk QCD medium. From a theoretical and experimental perspective 
heavy quark jets offer a rich new set of tools and challenges. Previous AdS/CFT calculations 
[10-19] of the heavy quark drag assumed a strongly coupled thermal medium of M = 4 super 
Yang-Mills (SYM) plasma. In this Letter we extend the earlier work on the heavy quark drag 
in a thermal medium to that of an arbitrary medium; in particular our formalism applies to 
both hot and cold strongly coupled nuclear matter. 

The AdS/CFT conjecture postulates a duality between certain field theories and the com- 
pactification of Type IIB string theory in various geometries [1-4] . In particular Af = 4 SYM 
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theory is dual to Type IIB string theory in AdS 5 xS 5 . What makes the conjecture so useful, but 
also so difficult to prove, is that the weak coupling limit of one theory is dual to the strong cou- 
pling limit of the other. Of especial interest to the heavy ion community is the strong coupling 
limit of QCD, where the only previous theoretical tool was numerical lattice simulations. While 
lattice simulations remain the only method to obtain quantitative results for strongly-coupled 
QCD, one may argue that for some observables similarities between QCD and M = 4 SYM 
theory can be exploited to improve our qualitative understanding of the former by performing 
calculations in the latter. In the limit of large 't Hooft coupling and number of colors N c , the 
M = 4 SYM theory is dual to an (often) analytically tractable theory, the classical limit of 
string theory: classical supergravity (SUGRA). 

Despite the many differences between QCD and M = 4 SYM, there have been a number of 
qualitative successes in applying the AdS/CFT ideas to heavy ion phenomenology at RHIC. 
Specifically, the lattice result of a surprising deviation of the entropy density above the decon- 
finement temperature from the Stefan-Boltzmann limit [20, 21] is easily understood analytically 
in AdS/CFT [22, 23]; the observed near perfect fluidity of the QGP [24-28] comes from a shear 
viscosity to entropy ratio an order of magnitude smaller than perturbative estimates, but in line 
with the AdS/CFT calculations [8]; the shockingly large suppression of heavy quarks [29-33] 
appears to be easier to describe in the AdS/CFT framework [10-13, 15-18]; and the emergence 
of an away-side double hump structure in two [34-36] and three particle [37, 38] correlators, 
while suggestive of a Mach cone formed by a supersonic jet [39, 40], can currently only be under- 
stood as the result of the near field energy and momentum deposition derived using AdS/CFT 
[19,41,42]. 

On the other hand, from both experimental and theoretical standpoints, perturbative meth- 
ods as applied to jet quenching physics should be viewed with increased skepticism. Early 
quantitative success in describing the energy loss of light partons [43-45] as evidenced by the 
suppression of pions [46, 47], etas [48], and the null control of direct gammas [49, 50] has 
given way to an inability [7, 51, 52] to simultaneously describe even any two of the following 
four observables: the (1) pion quenching and (2) similarly large suppression of non-photonic 
electrons (NPE), the decay fragments of heavy quarks, and the surprisingly large azimuthal 
anisotropy of the (3) hadrons from light parton jets [53, 54] and of the (4) NPE [31]. From a 
self-inconsistency standpoint the assumption of a small coupling, g <C 1, in high-p-r energy loss 
has never been true — at RHIC energies and temperatures g ~ 2 — with the momentum scale, 
set by either the Gyulassy-Wang model [55] of Yukawa-like scattering centers or the saturation 
scale, on the order of 0.5 — 1 GeV~ 0(1)Aqcd, far from perturbative. 

An understanding of the limits of applicability, of the inherent theoretical error, are required 
in order to either gain confidence in or falsify a calculation with experimental data. In deep 
inelastic scattering (DIS) and Drell-Yan production, the rigorous framework of factorization 
allows for a well-controlled pQCD expansion (see [56] and references therein). For heavy ion 
collisions it is not clear that factorization holds any longer [57]. Even supposing it does, the 
numerical value of the momentum characterizing the scale at which pQCD methods break down 
is not quantitatively known. The AdS calculations are under even less control. Until a dual 
string geometry is found for QCD one cannot quantitatively estimate the changes in observables 
from using a different theory; there is no path through dual theories approaching QCD whose 
parameterization may be used as an expansion parameter. 

What one can do is explore the results from different theories and geometries in an attempt 
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to discover universal behavior. This may give one reason to believe that the result will hold 
for QCD, at least at the qualitative level. In this Letter we follow this strategy by considering 
the heavy quark momentum loss in a novel geometry. Previous drag force calculations focused 
on motion of the gauge theory quark in vacuum (empty AdS metric) [15, 17, 18], or in a static 
thermal medium (black hole and black hole-like metrics) [10-13, 16, 17]. Here we examine the 
drag force on a static quark inside a shock wave. In particular the shock wave can be a model 
of the dense medium produced in a heavy ion collision. Equivalently it can also be a model of a 
nucleus incident on a probe quark, giving rise to cold-matter energy loss as is often considered 
in proton-nucleus collisions. 

The paper is structured as follows. In Sect. 2 we present the shock wave metric and solve 
the string equations of motion to find the trailing string solution, similar to [11—13]. We note 
that the analytic form of the trailing string solution, given below in Eq. (13), is considerably 
simpler than that found in [11-13]. In Sect. 3 we use the solution from Eq. (13) obtained in 
Sect. 2 to calculate the drag force on the quark inside the shock wave. Boosting the resulting 
expression to the rest frame of the medium and taking its typical momentum to be set by the 
temperature we obtain the same equation for the drag force as derived in the original works 
[11-13]. Just as the form of perturbative jet energy loss is independent of the state of the 
medium, it turns out that the heavy quark energy loss in a shock wave is exactly the same as 
that in a thermalized medium. Our result is derived for a shock wave made of any distribution 
of matter: due to Lorentz time-dilation the shock wave may represent a "snapshot" of any 
non-equilibrium matter. Hence our result generalizes the drag force of [11-13] to the case of a 
non-thermal medium. We also study the limitations of our calculation in Sect. 4 and find that 
the momentum "speed limit" of applicability for the shock wave calculation is parametrically 
exactly the same as for the black hole metric. We conclude in Sect. 5 by restating our main 
results and discussing the relation between the quark mass in vacuum and in medium for 
strongly-coupled gauge theories. 

2 Shock Metric, EOM, and the Trailing String Solution 

We will consider the generalized "shock" metric [9] 

L 2 

ds 2 = G^dx^dx" = ~2 [-2dx + dx~ + 2pz i d(x~) dx~ 2 + dx\ + dz 2 ] (la) 
L 2 

= — [- (1 -fiz 4 e{x~)) dt 2 -2pz 4 6(x')dtdx + (l + p z A 9(x~)) dx 2 + dx\ + dz 2 ] , 

(lb) 

where we have used the x ± = (t ± x)j \/2 normalization of light-cone coordinates with x = x 3 
and dropped the dVL\ standard metric of the five-sphere in AdS 5 x 5 5 . As usual L is the 
radius of the S 5 space. Also dx\ = (dx 1 ) 2 + (dx 2 ) 2 is the transverse part of the metric, and 
/i, v = 0, . . . , 4. 

The physical justification for describing Eq. (la) as a shock metric comes from the appli- 
cation of holographic renormalization [58], which relates the metric in Feffer man- Graham form 
[59] in d + 1 dimensions to the energy-momentum tensor in the d dimensional boundary theory. 
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For our particular choice of metric, Eq. (la), the above prescription yields 
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which describes an ultrarelativistic shock wave (shock front) moving in the positive x 3 = x 
direction. 

Previous calculations used /j,5(x~) as the coefficient of the dx~ 2 term in Eq. (la) to represent 
the Lorentz-contracted pancake of nuclear matter probed in DIS [60] or seen in heavy ion 
collisions [61]. Note that in those papers fi is a slightly different quantity with units GeV 3 , unlike 
GeV 4 here. As observed in [9] the coefficient of dx~ 2 can be any function of x~; the resulting 
metric still satisfies Einstein's equations. We take it oc 8(x~) to represent an incoming dense 
medium of nuclear matter colliding with a heavy quark in its rest frame. The setup is shown in 
Fig. 1, where the shaded region is the incoming shock wave. The fifth dimension z increases as 
one moves down the z axis: the top of the diagram is the boundary of AdS 5 space (z — 0); at the 
bottom [z = oo) is the stack of N c D3 color branes. A finite mass fundamental representation 
heavy quark in the 4D field theory corresponds in the 10 dimensional supergravity theory to 
the endpoint of an open Nambu-Goto string terminating on a D7 brane; the other end of the 
open string ends on the N c stack of D3 branes at z = oo. The D7 brane wraps an S 3 C S 5 
and fills the asymptotically AdS space from z = down to z = zm [62]. Before the collision 
the string hangs straight down (left panel of Fig. 1), while after the collision the string trails 
behind the quark (right panel of Fig. 1). 



Figure 1: Cartoon of the shock medium colliding with the heavy quark Q in its rest frame, (a) 
Before the collision the string representing the heavy quark hangs from the D7 brane straight down 
to the D3 branes at the origin, (b) After the collision the string trails behind Q. As the shock is 
moving in the positive direction (v s h > 0) in the heavy quark rest frame momentum actually flows up 
the string as the shock tries to transfer positive momentum to the quark; see Sect. 3 for more details. 

Ordinarily the spatial direction picked out when using light cone coordinates represents the 
beam direction. That would be the case here if we apply our model to the description of energy 
loss in proton-nucleus collisions, i.e., in cold nuclear matter. For jet energy loss in heavy ion 
collisions x corresponds to the direction of motion of the heavy quark in the lab frame, an 
orientation often taken transverse to the beam. 

While [1 1] extended the heavy quark drag calculations to more general metrics of the black 
hole type, a novel feature of the used here is its lack of an event horizon. If light can 
pass through a surface coming from inside the suspected black hole, the surface is not a true 
horizon. For our metric a light ray moving towards the boundary of AdS space both in the z 
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and x direction can cross the suspected horizon at Zh = /i -1 ^ 4 . Since the light ray can therefore 
escape from inside the suspected black hole, the surface at Zh = /i~ 1//4 is not a true horizon. 

We are interested in the motion of a heavy quark in the background specified by the metric 
of Eq. (lb). The test string action is 



Sng = -T J drday^g, 

g = detg ab , g ab = G^ d a X^ d b X u , 

where G^ u is the spacetime metric of Eq. (lb), Greek indices refer to spacetime coordinates, 
and Latin indices to worldsheet coordinates. X M = X^(a) specifies the mapping from the string 
worldsheet coordinates o a to spacetime coordinates x M . The backreaction of the fundamental 
string (0(N C )) is neglected as compared to the O(N^) contributions from the adjoint fields of 
TV = 4 SYM [63]. 

Varying the action, Eq. (3), yields the equations of motion: 

V a P; = 0, P^ = i^/ y /=g = -T G, a/ d a X v , (4) 

where V a is the covariant derivative with respect to the induced metric, g ab , and the 7r^ are the 
canonical momenta, 

d^g 

** = rn^r {) 

We limit our attention to only the nontrivial directions of t, z, and x, where x is parallel to 
the direction of the heavy quark propagation in the lab frame (although we will carry out our 
derivations in the rest frame of the heavy quark). X^(a) maps into the (t,x,z) coordinates; 
choosing the static gauge, a a = (t,z), the string embedding is described by a single function, 

x(t, z): 

X» = (t,0,0,x(t,z),z). (6) 

The equations of motion for x(t, z) can be derived in two equivalent ways: plug in the x(t, z) 
ansatz into Eq. (4); or first substitute the ansatz into Eq. (3), then vary the action. Following 
the latter approach and denoting d t x = x and d z x = x' we find 

-9 = ^(1 + x' 2 - x 2 - iiz\\ - xf), (7) 



d f [iz 4 — (1 + iiz^)x\ d 



x 



+ - — — =0. 



dt \ z 4 y^g J dz yz^y^g / 

One may determine the mass of the equivalent point particle in the 4D field theory by 
examining the string motion outside the shock. For \x = 0, the constant velocity string x(t, z) = 
xo + vt is clearly a solution to the equations of motion, Eq. (8); the action, Eq. (3), becomes 

S = -ToL 2 f dz\ f dtVT^P, (9) 

Jzm Z 2 J 

the action of a point particle of mass M q = T L 2 /zm- Hence we have that 



ToL 2 x/A 

ZM = -ALT = 2^' (10) 
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Figure 2: The three solutions, x(t, z) = £(z) = xq, xq ± y^Zz 3 /3, to the equation of motion for a 
static string inside the shock, Eq. (11). The constant solution (x(t, z) = xq) is unstable, S ~ ioo. We 
also throw away the unphysical time-reversed solution (x(t,z) = xq — ^//Iz 3 /3) by noting that any 
fluctuation would cause it to move toward the remaining (physical) one (see the text for details). 



where we have used the AdS dictionary to relate TqL 2 = y/\/2TV. 

Assuming an asymptotic static solution, x(t,z) = £{z), in a metric for which the shock 
occupies all space and time, the induced metric becomes time-independent. The equations of 
motion, Eq. (8), reduce to 



d ( £' 



dz \z 4 



0. (11) 



Consider the first integral of the motion resulting from Eq. (11). If we set the quantity in 
parentheses in Eq. (11) equal to a constant C, then we can solve for 



« z) = ±fe2 Vf=w (12) 

There are two cases to consider for a string hanging from the boundary down to the stack of 
D3 branes: C = and C ^ 0. (For other solutions of the equations of motion see Appendix 
A.) For C ^ 0, the constant of integration is fixed by considering the signs of the numerator 
and denominator inside the radical as a function of z: for small z both are positive; for large 
z both are negative. In order to avoid imaginary solutions, the numerator and denominator 
must change signs at precisely the same value of z; thus C = ^JJi. This leads immediately to 

£(z) =x ± ^z 3 , (13) 

which are shown in Fig. 2. It is amusing to note that the near-boundary expansion of the static 
quark solution for the black hole metric (with horizon at z = Zh) is x(t, z) m xq + vt ± z 3 / (3z^). 
However we would like to point out that the string configuration in Eq. (13) is valid for all z in 
the shock metric and is thus significantly simpler than that found in the black hole metric. As 
we will see below, while being analytically simpler, our solution retains all the main qualitative 
features of the black hole case. 

The sign ambiguity, resulting from the time-reversal symmetry of the problem, can be fixed 
by taking the positive, trailing solution; the negative sign has the string "trailing" in front 
of the heavy quark. We find from Sect. 3 that, in the lab frame, the physical solution loses 
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momentum to the medium whereas the unphysical solution gains it. Additionally we will find 
in Sect. 4 that any fluctuation of the negative sign string solution will lead to it approaching 
its time-reversed, physical solution. 

Since the heavy quark is static inserting Eq. (13) back into the action and integrating out 
the z coordinate yields a point particle action S = —M q J dt with the vacuum heavy quark 
mass M q . This should be contrasted to the result found in the black hole metric for which the 
point particle mass is modified from that in vacuum — interestingly, it is decreased in medium. 
We will discuss this further in the conclusions, Sect. 5. 

For C = 0, £ = xo is immediately found: the string hangs straight down. However plugging 
the solution back into the action gives 




One sees that the IR (large-z) part of the z integration gives the action an infinite imaginary 
part. We interpret this as an infinitely unstable state that would immediately decay into the 
physical trailing string solution. This and the two solutions from Eq. (13) are shown in Fig. 2. 



3 Drag Force 

The drag force on the heavy quark in the SYM theory corresponds to the momentum flow 
from the direction of heavy quark propagation down the string, i.e., dp/dt = —it]., it® are the 
canonical momenta and may be found from Eqs. (5) and (7): 

7T° \ T/ -4 / -l-X ,2 + /iZ 4 (l-x) 

< I =-^7= I ±-Hz\l-x) I (15) 



7T° J Z V 9 \y —x'[x — llZ A {l — x)) 



The "momentum change" of our heavy quark solution given by Eq. (13), where momentum 
change is in quotation marks as the quark is held static, is then 



dp i VA 

dt = -*'= 2^- W 

Since all our calculations have been in the heavy quark rest frame, Eq. (17) does not describe 
the momentum loss for a heavy quark in the lab frame, the rest frame of the shock. Even though 
formally the shock propagates at the speed of light we will think of it as an approximation 
to the physical setup of a medium moving almost on the light-cone; then this rest frame is 
well defined. Moreover we can now also relate /x to properties of the medium by comparing 
the energy-momentum tensor found from the metric using the holographic renormalization 
procedure, Eq. (2), to one derived from the medium properties. 

Following [60] , we assume the medium is made up of valence gluons of the M = 4 SYM 
fields; see Fig. 3. If in the rest frame of the medium the particles are isotropically distributed 
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with a typical momentum of order A — with associated inter-particle spacing of order 1 / A — then 
the 00 component of the stress-energy tensor in the rest frame of the shock is 

(Tu) cc iV c 2 A 4 , (18) 

where primes denote quantities in the rest frame of the medium and proportionality is up to 
a constant numerical factor. Changing into lightcone coordinates and boosting into the rest 
frame of the heavy quark yields 

(T__)=iV c 2 AV = iV c 2 A 4 (^) 2 , (19) 

where we assumed ultrarelativistic motion for the heavy quark in the medium rest frame, 
p' ~ Mj. (Eq. (19) can be viewed as a definition of the scale A absorbing numerical constants 

arising from the precise definition of (Tq ) and from the change to (T ).) Comparing this with 

Eq. (2) we read off 



H = 7T 2 A 4 



M 



(20) 



We now have the right hand side of Eq. (17) in terms of the momentum of the heavy quark in 
the rest frame of the medium; however the left hand side is still evaluated in the heavy quark 
rest frame. 




Vq = 




Vq = "V sh 
< 



Figure 3: Illustration of the shock wave colliding with a heavy quark in the (a) heavy quark rest 
frame and (b) medium rest frame. In (a) the shock moves with velocity v s h in the positive x direction; 
in (b) the heavy quark has momentum p' ~ jMq in the negative x direction. 



To rewrite dp/dt on the left of Eq. (17) in terms of the momentum and time in the medium 
rest frame, dp' /dt', note that dp/dt is the 3-vector component of the force 4-vector in the quark 
rest frame: 

r = ? = (2i) 



dr dt 

One can see from Eq. (16) that ~k\ = 0, and hence /* = 
shock rest frame is therefore 



-if 
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dt 



0, and the 4-force boosted into the 



(22) 



where the negative sign comes from boosting into a frame moving in the opposite direction; see 
Fig. 3. From the definition of the 4-force we also know that in this frame 



dp' 
dr 
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dp' 
dt'' 



(23) 
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Hence we find that 

dp = _<¥. (24 ) 
dt dt' 1 ; 

Eqs. (17), (20), and (24) lead us to the main result of this Letter: 

dt' ~ 2 M q P - [Zb) 

Should we take the typical medium particle momentum to be A = \^rT, reminiscent of the 
thermal equilibrium result, then our result exactly reproduces that of the black hole metric, 
dp'/dt' = -7iVATV/(2M g ) [11-13]. 

4 Limits of the Calculation 

As shown in [12, 13, 63] there are limits to the applicability of the heavy quark drag calculations 
in a black hole metric. In the original analytic construction the heavy quark was taken to have 
constant velocity, despite the momentum lost to the medium. If one proposes the existence of 
a field on the D7 brane supporting the finite mass heavy quark motion, then one is naturally 
led to a speed limit based on the maximum field strength supported by the Born-Infeld action. 
However numerical calculations were made for quarks with non-constant velocity and it turned 
out that the drag form derived analytically is a very good approximation for this more difficult 
problem. Nevertheless one might still have a speed limit induced by the warping of spacetime 
due to the presence of the thermal black hole. The finite heavy quark mass means that the 
string endpoint representing the quark in 4D propagates at z — zm in the 5 th dimension (see 
Eq. (10)). While at the Minkowski boundary, z = 0, the flatness of spacetime gives a local speed 
of light c = 1, reality of the point particle action in 5 bulk dimensions in the AdS BH metric with 
the horizon at z = Zh requires that \f—G^ v (dx^/dr) (dx u /dr) = yjl — (z / Zh) 4 ) / z 2 — v 2 /z 2 be 
real. Then v 2 < 1 — (z/z h ) 4 leading to 7 < [z h j ' z) 2 < (z h /z M ) 2 [12, 13, 63]. The speed limits 
resulting from these two lines of reasoning are precisely the same; this is not an accident but 
due to T-duality [64]. 

While the metric, Eq. (lb), does not support an event horizon, reality of the point particle 
action still yields an asymmetric, ^-dependent result for the local speed of light; to wit, in the 
rest frame of the heavy quark, 

<v<l. (26) 

At the boundary, z = 0, the usual Minkowski speed of light limit, —1 < v < 1, is recovered. 
Motion at the stack of D3 color branes (at z = 00) is restricted to the speed of light in the 
direction of the shock medium motion. 

In a similar way one finds that reality of the Nambu-Goto action requires the velocity of 
the string at z = 00 to be +1. This provides an argument against the reversed trailing string 
solution (given by Eq. (13) with the minus sign on its right hand side): the only classical 
fluctuations with non-zero velocity supported by it are those that, at the D3 branes at z = 00, 
give motion at the speed of light in the direction of the shock motion, i.e., towards the physical 
trailing string solution (given by Eq. (13) with the plus sign on the right hand side). This 
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demonstrates that the reversed trailing string solution is unstable under a large set of string 
fluctuations, and can therefore be discarded. (This conclusion should be compared to the 
stability analysis carried out in [65] for the black hole metric.) 

Plugging v = for our static quark into Eq. (26) gives us the bound for this calculation, 
namely 

fiz A M <l. (27) 

Using Eq. (10) to write z M = VX/(2nM q ) and Eq. (20) for p = vr 2 AV, along with A = y/nT, 
we obtain 



4ttM 2 4M 9 2 
A A 2 = AT 2 

The speed limit in this geometry is therefore identical to that for the BH metric [13, 63, 65] 



7 < q - = K (28) 

1 ~ AA 2 AT 2 1 ' 



5 Conclusions and Discussions 

In this Letter we used the AdS/CFT correspondence to calculate the drag force on a heavy 
quark in a strongly coupled medium. By placing the string in a shock metric, we extended 
the derivation of drag in a thermalized medium — whose metric contains a black hole — to a 
generalized medium not necessarily thermalized. It turns out that, just as for perturbative 
energy loss, the form of the drag is independent of whether the medium is thermal or not. 
Moreover the momentum speed limit for the applicability of our calculations in the shock 
metric is the same as was found for the black hole metric. 

To complete the comparison of our calculation with the existing ones for the drag force 
in a thermal medium [12, 13, 63] let us find the mass of the heavy quark in the shock wave. 
According to the standard prescription [12, 13, 63] the (static) quark mass in the shock is 

oo 

Mf ock = - J dzir° (29) 

ZM 

which, with the help of Eq. (15) and Eq. (13) yields 

ZM 

We see that the quark mass does not change inside the shock wave as compared to the mass in 
the empty space outside of the shock wave. This is in stark contrast to the thermal modification 
of the quark mass observed in [63] for a quark in a finite-temperature medium. 

To understand this result let us first discuss the origin of the medium modification of the 
quark mass found in [63]. * In [63] it was found that at finite temperature T the quark mass 

*We thank Cyrille Marquet for a discussion which led us to the qualitative understanding of the medium 
modification mechanism described below. 
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becomes M q {T) = (\/X/2tc) \{1/z m ) — (l/z h )] with the horizon radius z h = 1/nT. The mass 
of the heavy quark in AdS/CFT is constructed by integration over the length of the string, as 
shown in Eq. (29). On the gauge theory side this implies that the mass is built up from radiative 
corrections: at strong coupling the corrections are large and lead to a large dressed quark mass 
[66]. That is why the quarks in AdS/CFT are always heavy. Thinking of a heavy quark as 
being dressed by a cloud of strongly coupled particles ( "partons" ) it is easy to understand the 
field theoretic origin of medium modification of the quark's mass, due on the AdS side to the 
Zh IR cutoff in the ^-integration introduced by the BH horizon [17, 67]: the medium screens 
the infrared part of the parton cloud surrounding the quark, thus reducing the dressed heavy 
quark mass. It is essential that the quark is heavy not due to a large bare quark mass, but 
due to quantum dressing by the gauge field fluctuations: this allows the thermal medium to 
eliminate the soft IR fluctuations and make the quark mass smaller in the medium. 

Now, if we think of the shock wave considered here as of a boosted thermal medium (which 
is one of the possibilities), then the time scale of the thermal motion of the medium particles 
would be Lorentz-dilated in the rest frame of the quark: r t h ~ (V^) 7- As 7 ^> 1, the dilated 
thermal time scale is much longer than any other time scale in the problem. This is a 
well-known feature of the parton model in the infinite momentum frame, in which the partons 
appear "frozen in time" to the probe interacting with them. Hence our static heavy quark does 
not "feel" the thermal motion of the particles in the shock wave over the typical interaction 
time scale between the quark and a shock wave particle. This is why the heavy quark mass 
does not get modified by the shock wave, as we saw in Eq. (30): the lack of thermal motion 
on the interaction time scale means no screening for the infrared partons in the quark wave 
function. Indeed there is an IR screening due to saturation effects in the shock [60, 66], which 
may lead to quark mass modification if we integrate over z in Eq. (30) only up to the horizon of 
the induced metric on the string, as was done in [66, 67]. Our qualitative picture also explains 
why the thermal effects are not relevant for the quark energy loss: in the infinite momentum 
frame discussed here the thermal motion is irrelevant, but the quark still experiences a drag 
force. We conclude that the drag force indeed does not depend on the thermal motion, and is 
thus independent of whether the medium is thermal or not. This is exactly the same conclusion 
as we obtained by explicit derivation above. 
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A Other Solutions to the EOM 

The argument used to derive the stable static solutions x(t,z) = x , x ± ^/Jlz 3 /3 required the 
string to hang from a z smaller than both ^/Jl and \[C down to the stack of D3 branes at 
z = 00 so that, in order for the solution to remain real, both the numerator and denominator 
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of the radical in Eq. (12) had to change sign at precisely the same value of z. However strings 
with a turning point may also satisfy this reality condition. 

For such a solution x(z) to exist, z(x) should increase with x until the turning point, where 
the slope of z(x) must be (note the use of the inverse function). If we denote the value of z 
for which the numerator of the radical in Eq. (12) is zero as z c and the value of z for which the 



denominator is zero as z max , then 



4 = 1 - C'zt 

-1/4 _ (J-l/2 



V z t — 1 — C 2 z^ ax — (AT 



j max ■ 



=> z c = /J, 

The condition for turning, dz/dx = 0, thus requires the string solution to extend up to z = z n 

Concentrating on stable string solutions, ones for which the string configuration is purely 
real, we can consider two cases: (1) z M < z max < z c and (2) z M > z max , z c . The latter case 
may be solved analytically in closed form for the two shape possibilities: (a) string endpoints 
at plus and minus infinity in the x 3 direction, both terminating on the D7 or the stack of D3 
branes with a hump at xq that passes through z max and (b) a hanging string solution from the 
D7 brane to the D3 branes. There is considerable freedom in these solutions as C is no longer 
fixed. However none are physical due to the speed limit condition, Eq. (27). 

Unlike in the black hole metric one may also solve the hanging "meson" case (1) analytically. 
For simplicity let us consider infinitely massive quarks with zm = 0. The parameterization from 
x = —r/2 to for a string with ends at (x,z) = (±r/2, 0) may be found by direct integration 
of the first integral of the motion, Eq. (12): 



€(z) = -- + / dzz 




2 z 2 / \/ 1 — 

r z 3 ^ {3 1 1 7 4 z 4 . 
1 5 — Fi -; — ,-;-;uz 4 ,- — , (A.2 



max 



where Fx is Appell's hypergeometric function of the first kind (to construct the solution in the 
interval from x = to r/2 simply reflect Eq. (A.2) across x = 0); this solution is shown in 
Fig. 4. The solution in Eq. (A.2) complements the results of [60] in which the string profile was 
found for a dipole oriented perpendicular to the direction of shock wave motion, by providing 
the solution for a dipole oriented parallel to the shock velocity. 

Imposing the turning condition on Eq. (A.2), i.e., that in the middle of the string at x = 
one has the turning point z = z max , yields the following condition for z max : 

r v^Fr(7/4) _ f3 1 5 4 



2 3T(5/4) ZmaxF {r 2 ] A^ Zmax )- (A ' 3) 

Numerical analysis shows that for r > r crit w 0.85 z c = 0.85 /i -1 / 4 the solution for z max becomes 
complex- valued. 

For r < r cr it the singlet configuration for the quark anti-quark pair, Eq. (A.2), suffers 
no energy loss: the momentum lost by the leading heavy quark is exactly canceled by the 
momentum gain of the heavy quark that follows behind. However, for separations r > r crit the 
system is unstable to decay into a color-adjoint state, possibly leading to complete dissolution 
of the meson. To completely answer this question a full study of the roots for z max in the case 
of a general orientation of the dipole with respect to the shock velocity similar to what was 
done in [68] needs to be carried out: such analysis goes beyond the scope of this work. 
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Figure 4: Visualization of the hanging meson configuration, Eq. (A. 2). Note that z c > z max , and 
the string turns precisely at z = z max . 
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